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Almost universal codes achieving ergodic MIMO 
capacity within a constant gap 

Laura Luzzi and Roope Vehkalahti 


Abstract —This work addresses the question of achieving ca¬ 
pacity with lattice codes in multi-antenna block fading channels 
when the number of fading blocks tends to infinity. 

A design criterion based on the normalized minimum determi¬ 
nant is proposed for division algebra multiblock space-time codes 
over fading channels; this plays a similar role to the Hermite 
invariant for Gaussian channels. 

It is shown that this criterion is sufficient to guarantee trans¬ 
mission rates within a constant gap from capacity both for slow 
fading channels and ergodic fading channels. This performance 
is achieved both under maximum likelihood decoding and naive 
lattice decoding. In the case of independent identically distributed 
Rayleigh fading, it is also shown that the error probability 
vanishes exponentially fast. 

In contrast to the standard approach in the literature which 
employs random lattice ensembles, the existence results in this 
paper are derived from number theory. First the gap to capacity 
is shown to depend on the discriminant of the chosen division 
algebra; then class field theory is applied to build families 
of algebras with small discriminants. The key element in the 
construction is the choice of a sequence of division algebras whose 
centers are number fields with small root discriminants. 

Index Terms —MIMO, block fading, space-time codes, number 
theory, division algebras 


I. Introduction 

It is well-known (Tl] that in ergodic multiple-input multiple- 
output (MIMO) fading channels with channel state information 
at receiver only, the maximal mutual information is achieved 
with Gaussian circularly symmetric random inputs. In this case 
the existence of capacity-achieving codes can be proven with 
standard random coding arguments. 

It has been shown that by combining simple modulation and 
strong outer codes such as turbo or LDPC codes, it is possible 
to operate at rates close to capacity with small error probability 
|4j,l5t]. However, to the best of our knowledge, the problem of 
achieving capacity with explicit codes for all ranges of signal- 
to-noise ratio (SNR) is still open. 

This is in strong contrast to the classical complex Gaussian 
single antenna channel, where the capacity is log(l + SNR) 
and it is known that several lattice code constructions achieve 
log SNR— C rates. These constructions are based on a rich 
theory of lattice codes developed to attack these questions. 
At the heart of this theory are sphere packing arguments 
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that prove that the performance of a lattice code in the 
classical Gaussian channel can be roughly estimated by the 
size of a geometrical invariant of the lattice, the Hermite 
invariant. In particular the Hermite invariant can be used to 
roughly measure how close to capacity a family of lattices 
can get. This connection has been extremely fruitful and has 
led to a monumental work connecting algebra, geometry and 
information theory |@]. 

In the case of fading channels the situation is quite different. 
While it is well-known that space-time lattice codes from 
division algebras J3l provide good performance over multiple 
antenna fading channels, and a rich algebraic theory has been 
developed to optimize single codes [8|], there are as yet no 
results connecting capacity questions and the geometry of 
lattices. The minimum determinant criterion |9[] allows to 
improve the worst-case pairwise error probability in the high- 
SNR regime, when coding over a single fading block. Opti¬ 
mizing this value has been the major concern of several works 
in space-time coding iEI [Till . However, no design criterion 


has been suggested for approaching the MIMO capacity with 
explicit lattice codes. 

In this paper we address this problem and show that when 
we are allowed to encode and decode over a growing number 
of fading blocks, the normalized minimum determinant plays 
a similar role to the Hermite constant in Gaussian channels. 
In particular it can be used to measure how close to capacity 
a given family of lattice codes can get. 

Based on this design criterion we prove that for a MIMO 
channel with n transmit and n r receive antennas, there exists 
a family of multiblock lattice codes L n ^ C M nxn / C (C) 
(where k goes to infinity) that achieves a constant gap to 
capacity both in the slow fading and ergodic fading case. 
More precisely, for a MIMO channel with ergodic capacity 
C = E ff [logdet(7„ r + ^iytF)], our scheme achieves 
any rate 


R < E/y 


log det ^5; H^H 


n 


- n log C L +nlog—, (1) 

4 n 


where Cl is a certain geometric invariant of the family of 
lattices. These rates are achieved not only with maximum 
likelihood (ML) decoding, but also with naive lattice decoding. 
Furthermore, the same scheme achieves positive rates of reli¬ 
able communication for more general fading processes {Hi] 
under the mild hypothesis that the weak law of large numbers 
holds for the sequence of random variables {log det H^Hi}. 
As far as we know, this is the first explicit coding scheme 
which achieves constant gap to capacity for all SNR levels in 
MIMO channels. 
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Instead of using random coding arguments we consider 
algebraic multi-block division algebra codes [12, [13il . Our 
lattice constructions are based on two results from classical 
class field theory. First we choose the center K of the algebra 
from an ensemble of Hilbert class fields having small root 
discriminant and then we prove the existence of a K -central 
division algebra with small discriminant. Our lattices belong 
to a very general family of division algebra codes introduced 
in EHHQ. and developed further in lfl5|l and lfl6h . We 
will use the most general form presented in [17]. In particular 
our work proves that the classical number field codes HI 
achieve a constant gap to capacity in Rayleigh fast fading 
single antenna channels. 

In most works on algebraic space-time coding the code 
design criterion is derived from an upper bound for the 
pairwise error probability (s}] together with the union bound 
[II0. In our proofs we abandon this method and consider 
a hard sphere packing approach, classically used in lattice 
coding for the AWGN channel. The idea, formalized in Section 
EH is to exploit the special multiplicative structure of algebraic 
codes. It was observed in the context of diversity-multiplexing 
gain trade-off (DMT) analysis ifloi 21], that fading has a 


diminishing effect on the euclidean distance of the received 
code constellations derived from division algebra codes, hav¬ 
ing the so-called non-vanishing determinant property, only if 
the channel itself is bad. This property was formalized in 
[21], where the authors introduced approximately universal 
codes for fading channels. Our main results Theorem 14.11 
and Theorem 14.71 rely on this “incompressibility” property of 
algebraic lattices. It follows that our codes are almost universal 
and perform within a constant gap to capacity for a wide class 
of channels, having only mild restrictions on fading. 

While we discuss specific lattice codes from division alge¬ 
bras, our proofs do work for any ensemble of matrix lattices 
with asymptotically good normalized minimum determinant. 
The larger this value is, the smaller the gap to the capacity. 


This work also suggests that capacity questions in fading 
channels are naturally linked to problems in the mathematical 
research area of geometry of numbers. Unlike the single 
antenna Gaussian case, many of the questions that arise have 
not been actively studied by the mathematical community. 
Hopefully, studying such questions may lead to a comprehen¬ 
sive geometric theory of lattices for multiple antenna fading 
channels. 


We note that the proposed lattice code constructions are not 
yet practical, since they are based on number fields whose 
existence is proved through class field theory. Given a fixed 
degree, the required number fields can be found using compu¬ 
tational algebra software, but this process is computationally 
taxing. Decoding of the proposed codes is also very complex 
and the constructions we provide still have a large gap to 
capacity. 

On the other side, as demonstrated in Section IVIIII the 
existence results we use are very pessimistic. For small degrees 
the normalized minimum determinants of the best possible 
lattices are considerably better than the bounds provided by 
our existence results. 


A. Related work 


While our work shows that one can achieve a constant gap 
to capacity in ergodic MIMO channels with a fixed family 
of algebraic codes, it is natural to consider the more general 
question of whether it is possible to achieve capacity with 
any lattice codes. Such a result would be a generalization 
of the work in |22, [HI 24] which proved the existence of 


random lattice code ensembles achieving rate log(SNR) over 
the AWGN channel. By making the extra assumption that the 
transmitter and receiver have access to a common source of 
randomness in the form of a dither, the authors in [25|] finally 
proved that the AWGN capacity is achievable with random 
lattice codes. An explicit multilevel construction from polar 
codes was recently proposed in 126], 

As far as we know our work |2[] was the first to give a 
proof that lattice codes achieve a constant gap to capacity in 
block fading MIMO channels. In the single antenna fast fading 
channel this problem was considered before in (27^, which 
claims that random lattices achieve a constant gap to capacity. 
In |[28| the authors extend their previous results and claim to 
give a proof that random lattices achieve capacity in single 
antenna ergodic fading channels. However, we believe that at 
least in its current form, the analysis in both works is missing 
some fundamental details. In particular, the gap A < 1 + 
logE/j tttj , given in lEH Theorem 3], is infinite even when 


the fading process {hi} is i.i.d. complex Gaussian. In 
Equation (20)] the authors state that for a given fixed fading 
realization, the Minkowski-Hlawka theorem implies that there 
exist a lattice for which the error probability is upper bounded 
in a certain way. However, they proceed as if there existed 
a single lattice that would satisfy this upper bound for any 
channel state. To the best of our knowledge, such a result can 
not be derived from Minkowski-Hlawka. 

As far as explicit algebraic constructions are concerned, our 
work is indebted to several previous papers. 

The idea to use division algebra codes to achieve capacity 
can be tracked down to the work of H.-f. Lu in 1131] . While 
studying the diversity-multiplexing gain tradeoff (DMT) of 
multiblock codes he conjectured that the ensemble of multi¬ 
block division algebra codes might approach the ergodic 
Rayleigh fading capacity. Our work confirms that conjecture; 
however, we point out that it is unlikely that DMT-optimality 
alone is enough to approach capacity. Instead one should 
pick the code very carefully by maximizing the normalized 
minimum determinant. 

The families of number fields on which our constructions 
are based were first brought to coding theory in IHI , where 
the authors pointed out that the corresponding lattices have 
large Hermite constant. C. Xing in 130] remarked that these 
families of number fields provide the best known normalized 
product distance making them a natural candidate for achiev¬ 
ing constant gap to capacity in fading single antenna channels. 

Our results on slow fading channel are motivated by the 
work in 1311] . where the authors prove that precoded integer 
forcing achieves a constant gap to capacity for every slow 
fading channel with fixed fading. 


Our geometry of numbers approach has its roots in 13211 , 
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where the authors studied lattice codes in single antenna fading 
channel and defined the normalized product distance. They 
also pointed out that using this criterion reduces the lattice 
design to a problem in geometry of numbers. 

The generalization of these ideas to the MIMO channel 
was developed in [33] and 1 13411 . where the code design for 
quasi-static MIMO channel problem was translated into lattice 
theoretic language and where a formal definition of normalized 
minimum determinant was given. However, none of these 
works considered the relation between geometry of numbers 
and capacity problems. 


B. Organization of the paper 

In Section QT| we introduce the multiblock channel model 
and recall the relevant properties of lattice codes. In Section 
IHI-AI we develop a geometric design criterion for capacity 
approaching lattice codes for fading multiple antenna channels 
and define the concept of reduced Hermite invariant which is 
an analogue of the classical Hermite invariant. In Section ITlI-B I 
we state the existence of lattices having asymptotically good 
normalized minimum determinant (the proof will be given in 
Section [Yin . In Section [IV] we prove that the lattice codes of 
the previous section achieve positive rates over a very general 
class of channels. We then prove that they achieve a constant 
gap to capacity over slow fading (Section IV-Ab and ergodic 
fading channels (Section [V-BI >. In Section[VI] we focus on the 
i.i.d. Rayleigh fading channel model, and show that the error 
probability vanishes exponentially. In Section lVTIl we prove the 
existence of asymptotically good lattices, and in Section lYIHI 
we specialize our results to the single antenna case. Finally 
in Section [IX] we explore the connection between capacity 
questions in fading channels and geometry of numbers. Section 
[X] discusses some perspectives and open problems. 


C. Notation 

Throughout the paper, capacity is measured in bits. Ac¬ 
cordingly, we denote by log the base 2 logarithm in rate and 
capacity expressions; the natural logarithm will be denoted by 
In. 


II. Multiblock lattice codes 
A. Channel model 

We consider a MIMO system with n transmit and n r receive 
antennas, where transmission takes place over k quasi-static 
fading blocks of delay T = n. Each multi-block codeword 
X £ M nxn k( C) has the form (Xi, X 2 , ■ ■ ■, Xk), where the 
submatrix Xj £ M n (C) is sent during the i-th block. The 
received signals are given by 

y, = H,X l + W {, i£{l,...,k} (2) 

where Hi £ M„ rXn (C) and Wj £ M nrX T(C) are the channel 
and noise matrices. The coefficients of Wi are modeled as 
circular symmetric complex Gaussian with zero mean and unit 
variance per complex dimension. Perfect channel state infor¬ 
mation is available at the receiver but not at the transmitter, and 
decoding is performed after all k blocks have been received. 


We will call such a channel an ( n,n r: k)-multiblock channel. 
In this paper, we will also assume that for all i > 1, Hi £ 
M rir x n is full-rank with probability 1, and that the random 
variable i log det(HjHi) converges in probability to 

some constant when the number of blocks k tends to infinity. 
This channel model covers several standard MIMO channels 
such as the Rayleigh block fading channel and the slow fading 
channel. 

A multi-block code C in a (n,n r ,k )-channel is a set of 
matrices in M„ xn j.(C). In particular we will concentrate on 
finite codes that are drawn from lattices. Let R denote the 
code rate in bits per complex channel use; equivalently, |C| = 
2 Rkn . We assume that every matrix X in a finite code C C 
M nxn k(C) satisfies the average power constraint 

\ imi 2 < p, (3) 

nfe 

where ||X|| is the Frobenius norm of the matrix X. 


B. Lattice codes 

Definition 2.1: A matrix lattice L C M nxn k( C) has the 
form 

L = ZBi © ZB 2 © • • • 0 Z B r , 

where the matrices B\ 1 ... ,B r are linearly independent over 
R, i.e., form a lattice basis, and r is called the rank or the 
dimension of the lattice. 

The space M nxn k{ C) is a 2?r 2 fc-dimensional real vector 
space with a real inner product 

{X,Y)=K{Tr{XY')), 

where Tr is the matrix trace. This inner product also naturally 
defines a metric on the space M nx nfc(C) by setting \\X\\ = 

Given an m dimensional lattice L C M nXn k(C), its Gram 
matrix is defined as 


G(L) = ((X i ,X j )) 1 ^ m , 

where {3fi}i<i< m is a basis of L. The volume of the 
fundamental parallelotope of L is then defined as Vol(L) = 

V\ det(G(Z/))|. 

In the following we will use the notation R (L) for the linear 
space generated by the basis elements of the lattice L. 

Lemma 2.2: [35] Let us suppose that L is a lattice in 
■^nxfen(C) and S' is a Jordan measurable bounded subset of 
R(L). Then there exists X £ M nx /. n ( C) such that 


|(L + x)ns| > 


Vol(S) 
Vol(L) ‘ 


Given a family of lattices L n ^ Q M„ X nfc(C), let us now 
show how we can design multiblock codes C having rate 
greater or equal to a prescribed constant R. and satisfying 
the average power constraint (0, from a scaled version aL n ^ 
of the lattices, where a is a suitable energy normalization 
constant. We denote by B[r ) the set of matrices in M nxn k(C) 
with Frobenius norm smaller or equal to r. According to 
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Lemma [Z2l we can choose a constant shift Xr £ M nxnk ( C) 
such that for C = B(y/Pkn ) D (X# + aL n , k ) we have 

= |r| > Vol (BjVPkn)) = C n , k P n2k 
M - Vol (aL n , k ) Q,2n 2 fc Vol(L nj fc) ’ 


, nrfe 

where • We ^ en find ^e following condition 

for the scaling constant: 


a 


2 



Vol(L„, fc )^ 


(4) 


III. Design criteria for fading channels 

In this section we propose a new design criterion for 
capacity approaching lattice codes in fading channels. For the 
sake of simplicity, we will focus on the case n r > n. We 
note that the design criterion derived here will finally be the 
familiar minimum determinant criterion. However, we hope 
that our alternative characterization offers more insight on the 
topic and can have applications in further research. 


A. Reduced Hermite invariant 

We recall the classical definition of the Hermite invariant, 
which characterizes the density of a lattice packing: 

Definition 3.1: The Hermite invariant of an m-dimensional 
lattice L C M nxnk (C) can be defined as 

wm inf{||X|| 2 |XgL,X^0} 

( ) Vol(L) 2 / m 

On the n x n MIMO Gaussian channel such that the channel 
matrices Hi = I n Vi, the classical sphere packing approach is 
to choose a 2 ? 7 2 fc-dimensional lattice code L Utk C M nxnk (C) 
such that h(L ntk ) is as large as possible. 

Let us now assume that we have a finite code Cl C 
L nk C M nxnk (C) and a random channel realization 
H = [7/ |..... H k ], We note that the channel output Y = 
[Yi,..., Y k \ can be written as 

Y = H d X + W, 

where X = [X 1 ,...,X k ], W = [Wi,...,W fc ], and H d = 
diag(L/i,... ,H k ). From the receiver’s point of view, this is 
equivalent to an additive white Gaussian noise channel where 
the lattice code is 

HC l = {H d X | X e C L }. 

Even if the lattice L rt k (and therefore the code Cl) has good 
minimum distance, there is no guarantee that the same can be 
said about the lattice HCl■ This leads us to consider matrix 
lattices L nk C M nxnk ( C) which would have good minimum 
distance after any (reasonable) channel. If we assume that each 
of the matrices Hi in equation © has full rank with probability 
1, then the multiplication X K > H d X is a bijective linear 
mapping with probability 1. For any lattice L n k C M nxnk ( C) 
having basis Hi,..., B 2n 2 k we then have that 

HL n k = {H d X | X £ L Uik } = ZH d Bi ® ■ ■ • © 7iH d B 2n 2 k , 


is a lattice with basis H d B i, • • • , H d B 2n 2 k , and h (HL n , k ) is 

well defined. 

As a discrete group, HL n k has positive Hermite invari¬ 
ant, but even if h(L n ,fc) is large there is no guarantee that 
h( HL n k ) is. 

For convenience we first introduce a group of matrices 

G = {H £ M nXnk (C) | pdet(iT) = 1}, (5) 

where pdet (H) = ]^ =] det(Hi). 

Definition 3.2: The reduced Hermite invariant of an 771- 
dimensional lattice L C M nxnk (C) with respect to the group 
G is defined as 


rh G (L) = inf {h(HX)}. 

rl GCj 

For any lattice L, h (L) > 0. The same is not true for 
the reduced Hermite invariant. Let us now describe the set 
of lattices L for which rh G (L) > 0. 

Definition 3.3: The minimum determinant of the lattice L C 
M nxnk ( C) is defined as 


det m i n (L) '.— 


inf I pdet (X) I. 

A'eL\{0} 


If det m ,;„ ( L ) > 0 we say that the lattice satisfies the non¬ 
vanishing determinant (NVD) property. 

We can now define the normalized minimum determinant 
S(L), which is obtained by first scaling the lattice L to have 
a unit size fundamental parallelotope and then taking the 
minimum determinant of the resulting scaled lattice. A simple 
computation proves the following. 

Lemma 3.4: Let L be an m-dimensional matrix lattice in 
M nX nfc(C). We then have that 


St j \ tlet mi n (L) 

~ (Vol (L)) nk / m ' 


( 6 ) 


The normalized minimum determinant provides an alter¬ 
native characterization of the reduced Hermite invariant, but 
before that we need a well known lemma. 

Lemma 3.5: Let A be an m x m complex matrix. We have 
the inequality 


| det (A) | < 


m m / 2 


For a matrix X £ M nxnk ( C) this immediately implies that 

™ ^ 

Proposition 3.6: If L C M nxnk (C) is a 2n 2 /c-dimensional 
lattice, then 

nk (S(L)) 2/nk = rh G (L). 

Proof: If the lattice L includes a non-zero element X 
such that pdet(X) = 0, it is easy to see that nk (S(L)) 2 ^ nk = 
rh G (L) = 0. 

Let us now assume that pdet(X) 7 ^ 0, for all X ^ 0. If 
pdet ( H) = 1, Lemma [331 implies that 

||#dV || 2 > nk \pdet(H d X)\ 2 ^ nk = nk |pdet(X)| 2 ^ nfe . 

It follows that nk (8(L)) 2 ^ nk < rh G (L). 
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Let us now assume that we have a sequence of codewords 
iW g i suc h that 

lim nfc|pdet(JV^)| 2 ' / ” fc = nk (S(L)) 2 ^ nk . 

i—> oo 

If X W = [x[ l \ ..., X^], we can choose = 

pdet(XW) 1 / n [(X 1 (i) )- 1 ,..., (X^)- 1 ] so that pdet(ffW) = 
1. We then have 

\\Hfx ^\\ 2 = nk\pdet{Hf X^)\ 2 / nk = 

= nk\pdet(X^)\ 2 / nk 

for every i and therefore 

lim \\Hfx^\\ 2 = nk (S(L)) 2/nk . 

i—too 

It follows that nk (S(L)) 2 ^ nk = rhc(L). □ 

Remark 3.7: Our definition of the reduced Hermite invari¬ 
ant rh <3 depends heavily on the group G. The group chosen 
in © can be seen as a block diagonal subgroup of SLfc n (C). 
We could also consider a subgroup G\ C G and define rh ( ;; 
with respect to this group. A natural consequence of these 
definitions is that for two subgroups G i, G 2 of G such that 
G\ C G 2 , we have that 

rhGi(i) > rh G2 (L). 


B. Asymptotically good families of lattices 

Based on the observations in the previous section, we can 
say that a sequence of codes L n ^ is asymptotically good for 
the AWGN channel if lr(L ra ,fe) > cn 2 k, for some positive fixed 
constant c. Similarly, we can say that a sequence of lattices is 
asymptotically good for fading channels if rho (Tn,fc) > cn 2 k. 
As seen in Proposition 13.61 this is equivalent to asking that 
8{ L n,k) 2/nk > cn. 

In order to keep the paper suitable for a larger audience we 
will postpone the proof of the following existence result to 
Section IVnl 

Proposition 3.8: Given n, there exists a family of 2 n 2 k- 
dimensional lattices L n ^ C M nxn k(C), where k grows to 
infinity, and a constant G < 92.4 such that 

kn(n-l) ( (7\ ” 2fc 

Vol(L„, fc ) < 23 10 ( -J 

det min{L n ,k) = 1 and S(L nJ c ) > 23 ^ (n _ 1) ^ G y 2 ^ • 

Remark 3.9: In this section we have developed the notion 
of reduced Hermite invariant for the case n r > n. We observe 
that this notion does not extend in a straightforward way to the 
case n r < n, because the image HL of an infinite lattice L 
will no longer be a lattice, and the minimum distance in HL 
will be zero. However, when considering finite constellations 
Cl, it is still possible to find suitable lower bounds on the 
minimum distance of the received constellation HCl, as will 
be shown in the following sections. 


IV. Achievable rates for general channels 

Suppose that we have an infinite family of lattices L). £ C k 
with Hermite invariants satisfying > c, for some positive 
constant c. Then a classical result in information theory states 
that with this family of lattices, all rates satisfying 

( 7TC \ 

—J + logc, 

are achievable in the additive complex Gaussian channel fl, 
Chapter 3]. This means that we can attach a single number 
h(Lfc) to each lattice Lk £ C fc , which roughly describes 
its performance and in particular estimates how close to the 
capacity a family of lattices can get. The following theorem 
can be seen as an analogue of this result for fading channels. 

Theorem 4.1: Suppose that n r > n, and let {Hi}i e z be 
a fading process such that Hi £ M„ r . x „ is full-rank with 
probability 1, and suppose that the weak law of large numbers 
holds for the random variables {logdet(fT|fTi)}, i.e. 3p > 0 
such that Ve > 0, 


lim P 

k—> 00 


^J2logdet(HjHi) - p 



(7) 


Let L n & C M nx „fe(C) be a family of 2?r 2 /c-dimensional 
multiblock lattice codes such that 


det„ 


(L n , k ) = 1, and Vol (L n , k )^ < C L (8) 


for some constant Cl > 0. Then, any rate 

R< p + n (log P - log C L + log 

is achievable using the codes L n ^ both with ML decoding 
and naive lattice decoding. 

Remark 4.2: We note that existence of a family of lattices 
with 

C L < 23^ (j) , 

was given in Proposition 13.81 

Remark 4.3: This theorem is stated by giving two condi¬ 
tions © for the lattices However, according to Lemma 
13.41 we could have captured both of these conditions by an 
equivalent assumption 5(L n: k) 2 ^ nk > Proposition 13.61 
then transforms this condition to 

rho (L n ,k) = nk (5(L nt k)) 2/nk > 77 

L L 

As only k is growing, we can further write that rhG(Ln,fc) > 
n 2 kC' L , where C' L = n/CL■ We can therefore see that 
conditions © assure that the family of lattices L n j. is asymp¬ 
totically good in the sense of Section UlI-BI 

The achievable rate R in Theorem 14.11 can then be seen as 
a complete analogue to the classical sphere packing result in 
AWGN channels. 

Remark 4.4: The condition (© holds in particular for er- 
godic stationary fading channels and for constant MIMO 
channels. These special cases will be analyzed further in 
Section El where we will show that the codes in Theorem 
14.11 achieve a constant gap to channel capacity. 

To prove Theorem 14. II we need the following Lemma: 
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Lemma 4.5: Consider the finite code C = B(\/Pkn) D 
(Xr + aL n ,k) defined in Section III-BI Suppose that the 
receiver performs maximum likelihood decoding or “naive” 
lattice decoding (closest point search in the infinite lattice). 
Then, under the hypotheses of Theorem 14.11 Ve > 0 the error 
probability is bounded by 


P e < 2e~ 


a 
4 n 


JQ det(iT| Hi) nk < 1 + e l (9) 


P, < 


iw 

kn 2 


>i + e +: 


'{iS? <1 + e 


y kn 2 


^ 1 + 6 / < 2 6 


a 
4 n 


JJ \det(Hi)\ nk <1 + 6 


i= 1 


and Rj = 


R'i 

o 


Proof: We distinguish two cases: the symmetric case 
where n r = n, and the asymmetric case with n r > n. Let 
(In denote the minimum Euclidean distance in the received 
constellation: 


d H = min S \\Hi(Xi - A.;) 
x,xec -f—' " 

XjtX 2-1 


, with R[ £ M n ( C) upper triangular. Note 

that the “thin” QR decomposition Hi = Q'iR'i also holds. 
Multiplying the channel equation (Q) by Q], we obtain the 
equivalent system 

Yi = QjYi = RiXi + QjWi 

for all i = 1 ,k. Note that 


Yi = 


Thus, the second component contains only noise and no 
information. The output of the naive lattice decoder can be 
written as 


[ V. 1 


' R'iXi + Q\Wi ' 

L y c \ 


Q”Wi 


X = argmin £ \\Yi - H z X'f = 

X f £oiL r . k • i 

n ’ K i=l 


a) Case n r = n: We have 

Pe<Pl\\W\\ 2 >(^- 


where W = (W±, W 2 , ■ ■ ■, W i- ) is the multiblock noise. By 
the law of total probability, Ve > 0 we have 


( 10 ) 


= argmin Y) — RiX[ 

X l &aL nk Z ~f" 

n,/c l= i 

k 

= argmin £ (|| Y' - R^X'f + ||Q"W,|| : 
X'eaL n , k v 

k 

= argmin y \\Y( - R^X-f , 


Note that 2 ||kL|| 2 ~ x 2 {2kn 2 ). and the tail of the chi-square 
distribution is bounded as follows for e £ (0,1) 


( 11 ) 


Thus, the first term in equation ( ITOt vanishes exponentially 
fast as k —> 00 . 

In order to provide an upper bound for the second term, 
we consider a lower bound for the minimum distance in the 
received constellation. We have 

k 

d 2 H > a 2 nk min TT |det(fTiAj)| > 
xrL.Alo } 11 
2—1 

k 

> a 2 nfc |det(7Tj)|" 7? , 

i=i 

where the first bound comes from Lemma [331 and the second 
from the hypothesis that det m i n (L n ^) = 1. Therefore, the 
second term in (TITTb is upper bounded by 


X’£uL n 

since the second component does not depend on the lattice 
point X'. Thus, the naive lattice decoder for the original 
system declares an error if and only if the naive lattice decoder 
for the (n, n, k) multiblock system with components 

Y( = R’X t + Q\Wi = R’X, + W- 

does. Note that W[ = (QtyWi is an n x n matrix with i.i.d. 
Gaussian entries of variance 1 per complex dimension. 

Let dn' be the minimum distance in the 2n 2 fc-dimensional 
constellation generated by R' = \If ...., R' k \: 


d.R> = nnn Y R Xi-Xi 
x,x ec“ 

XrjLX i_1 

Observe that Vi = 1,..., k, 


( 12 ) 


| \Hi(Xi - A. t )|| 2 = || Q\K{Xi - Xi)|| 2 = || R'fXi - Xi )11 2 

Thus, dn = d,R. Moreover, det(+//,;) = det {{R'f)^R[) = 
|det(f?')| 2 . Similarly to the symmetric case, the error proba¬ 
bility can be bounded by 


P e <P( IIVL'II 2 > 1 dR ' 


b) Case n r > n: In this case, the lattice HdL n ^ is 2 n 2 k- 
dimensional but is contained in a 2n r n/.:-dimensional space. 
For all i = 1 ,,k, consider the QR decomposition 

Hi = QiRi, Qi £ (Hn r xn r (€)) Hi C. M nr xn(fd-) i 

where Qi is unitary and H, is upper triangular. We have Q, = 
[Qi Q"], where Q[ £ M nrXn (C) is such that {QtyQ'i = I n , 


where W' = , Wf). We can write 

k k 

d 2 R f > a 2 nk |det(f?')| nk = a 2 nk det(JT|fTj)^. 


i =1 


2=1 


The proof then follows exactly the same steps as in the 
symmetric case. □ 
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Proof of Theorem 1 4.1\ The second term in <[9]» can be 
rewritten as 


V. Achieving constant gap to capacity for slow 

FADING AND ERGODIC CHANNELS 



X logdet (HjH,) < log 




and will vanish as long as 


log 


( 4n(l + e)\ 

l « 2 ) 



A. Slow fading channel 

We now consider a slow fading scenario, where Hi = H 
is constant. When H is known both at the transmitter and 
receiver, the channel capacity is given by 

C(P) = max logdet(P + HQ X H^), 

Qx>0,tr(Q x )<P 


Recalling that the normalization constant a 2 in equation ((TJi 
satisfies 

r~il/n 2 k p 

Q 2 > 

- 2 R/*C l 

under the hypothesis that Vol(L n ,fe)"^ < Cl, a sufficient 
condition to have vanishing error probability is 

— < log P + - - log(4n(l + e)) - log C L + -2— log C n>k 
n n n z k 

From Stirling’s approximation, for large k we have 

(C n , k )Tk ^ne/{n(2irn 2 k)^). (13) 

Since log 27m 2 fc -A 0 when k -A oo, any rate 


R < /t+n(logP—log(4n(l+e))—logCL+log7re—logn) 


is achievable. This holds Ve > 0, and concludes the proof. □ 

Remark 4.6: Note that the two-sided convergence in prob¬ 
ability in equation 0 is actually not required in the proof of 
Theorem 14. II The theorem still holds provided that Ve > 0, 

lim P < (i — y y log det (Hj HA > el = 0. (14) 

^°° { ) 

Moreover, if we have exponentially fast convergence in (IT4l) . 
then the error probability P e also vanishes exponentially fast 
when k —> oo. 

As a final remark, we note that we can prove an analogue 
of Theorem 14.11 also in the case n r < n, although the bound 
on achievable rates is more involved: 

Theorem 4.7: Suppose that n r < n, and let {7p}j 6 z be 
a fading process such that H, £ M n , r x „ is full-rank with 
probability 1. Suppose that the weak law of large numbers 
holds for the random variables logdet (HiHj), i.e. 3p > 0 
such that Ve > 0, 


lim P 

k—too 


^ XIlogdet (HiH^)-n 



(15) 


Let L nj fc C M nxn fc(C) be a family of 2n 2 fc-dimensional 
multiblock lattice codes satisfying ®. Then, any rate 

7T6 

R < p + n r (log P — 2) + (n — n r ) login - n r ) + n log 

n z C L 

is achievable using the codes L n , k both with ML decoding 
and naive lattice decoding. 

The proof of Theorem 14.71 can be found in Appendix lAl 


where Q x is the covariance matrix of the input x for a single 
channel use. 

However, if the channel is known at the receiver but not at the 
transmitter, the transmitter cannot use optimal power allocation 
and waterfilling, and can only achieve the white-input capacity 
corresponding to uniform power allocation Q x = P/ n : 

Cwi = logdet ^ 'I nr + = logdet . 

This is for example the case for an open-loop broadcast 
channel where the transmitter cannot perform rate adaptation 
for all the users. 

Clearly, Theorems 14.11 and 14.71 apply to the slow-fading 
scenario since the law of large numbers holds. Moreover, 
the convergence of the error probability to zero will be 
exponential, since the second term in equation 0 is actually 
zero. The following corollary then shows that a constant gap 
to white-input capacity is achievable: 

Corollary 5.1: Consider a slow fading channel such that 
Hi = H for all i 7* 1. and let CZ be a 

family of 2n 2 fc-dimensional multiblock lattice codes such that 
det min (L n>fc ) = 1 and Vol(L njfe )^ < Cl- Then, this coding 
scheme can achieve any rate 

P 7T6 

R < log det —H^H — n log Cl + n log — 
n 4 n 

if n r > n, and any rate 

R < logdet —HH^ — 2n r —(n — n r ) log-hnlog 2, 

n n — n r uCl 

if n r < n. 

Remark 5.2: In the case n r > n, let A,;, i = 1 
be the singular values of H. Then the channel capacity can 
be written as C(P) = XlLi 1°§ (l + -^A*) . The previous 
corollary shows that the achievable rate is of the form 

R(P) = max [ 0, log det — jp^H — c 
V n 

for some constant c > 0. Let P m i n be the smallest value of P 
such that R(P) > 0 if P > P m ; n . Then, for P < P m i n we have 
that C{P) - R(P) = C(P) < C(P min ), while for P > P min , 
C(P) — R{P) = Y2i= 1 1°S (l + tP) + c which is a strictly 
decreasing function of P and tends to c when P —► oo. Thus, 
for all P > 0 we have that C(P) — R(P) < C(P m i„). This 
shows that the gap is bounded by a constant, however the value 
of P m in and of the constant depends on the channel. 

A similar argument holds for n r < n. 










B. Stationary ergodic channels 

We now specialize the results of Section [IV] to the case 
where the fading process {//,} is ergodic and stationary. For 
the sake of completeness, we review the relevant definitions 
here. 

Let I be the set Z or N, and consider a random process 
X 1 = {Xi} i 6 x on a probability space where each 

random variable X, takes values in a separable Banach space 
X. The sequence space {X 1 , B(X x )) with the Borel sigma- 
algebra inherits a probability measure mx from the underlying 
probability space, defined by 

mx(A) = P {w : X x (u>) e A} \/A £ B(X x ). (16) 


Definition 5.3: The process {Xi} is called stationary 
if Vi, k £ N, \/ii, * 2 , • • •, ik G I, the joint distri¬ 
bution of (Xjj, Xi 2 ,..., Xi k ) is the same as that of 

{Xi 1+t , x i2+t ,..., x ik+t ). 

In this case it is well-known [37], p. 494] that the measure 
mx is invariant with respect to the shift map T : X 1 —>■ X x 
such that T({xj}) = {tCi+i}. 

Definition 5.4: The process { X ,} is called ergodic if MA £ 
B{X x ) such that T~ 1 {A) = A, we have that mx(A) is equal 
to 0 or 1. 

We now go back to the channel model ({2}. For the sake of 
simplicity, we suppose that n r > n. If the fading process {//,} 
is stationary and ergodic, it is not hard to see that the random 
process {A'i} = |logdet(7/l//i) j taking values in X = R 
is also stationary and ergodic, and the shift T : R 1 — > R 1 
preserves the measure mx defined in ( IT6t . 

For an ergodic process such that the shift T is measure¬ 
preserving, Birkhoff’s theorem |j38 'l guarantees that for any 
f £ L 1 ^ 1 , B{X x ), mx), the sample means with respect to 
/ converge almost everywhere: for almost all {xj} £ X x , 


n=l 




fdm x - 


(17) 


In particular, the projection II : R 1 —> R on the first coordinate 
is L 1 according to the image measure mx if and only if 
E [|logdet7/t7/|] < oo. Under this hypothesis, Birkhoff’s 
theorem implies the law of large numbers: 

k 

lim ~y^Xi= [ I\.({xi})dmx({xi}) = 
fe_>0 ° k i=l 7h x 

= / noX I dP= I XidF = E[X] a.e. (18) 

J f2 J fl 

In other words, 

k 

lim J y2 log det(HjHi) = E H [log det(7/l//)] (19) 

k —y oo n " 


i= 1 


almost everywhere. 

In the ergodic stationary case, it is well-known [ 3, 39] that the 
ergodic capacity of the channel is well-defined and does not 
depend on the channel correlation with respect to time, but 
only on its first order statistics. Given a power constraint P in 


equation ([3]), the ergodic capacity (per channel use) is equal 
to 


C(P)= max E h [logdet(/„ r + HQ X H^)\ , 

Q x >0,tr(Q x )<P ^ J 


where H is a random matrix with the same first-order distri¬ 
bution of the process {Hi}, which is independent of time by 
stationarity, and Q x is the covariance matrix of the input x 
for one channel us^3. 

If we suppose that the channel is isotropically invariant , i.e. 
the distribution of H is invariant under right multiplication by 
unitary matrices, then the optimal input covariance matrix is 
<3x = £ In [§] and we have 


C{P) = Ep 



Since det(/ nr . + -^TTITl) = det(/ n + we can also 

write 


C(P) = E ff 




The following Corollary to Theorem 14.11 shows that in this 
case, the proposed multiblock codes can achieve a constant 
gap to ergodic capacity. 

Corollary 5.5: Suppose that n r > n and that the fading 
process {Hi} is ergodic, stationary and isotropically invariant. 
Moreover, suppose that E [|logdet < oo. Let L n ^ C 

Af„ X nfc(C) be a family of 2n 2 /c-dimensional multiblock lattice 
codes such that det m i n (L n ,fc) = 1 and Vol(L n ,fe) ^ Cx. 
Then, any rate 


R < E h 


log det —H^H 
n 


- n log Cl + n log 


7re 
4 n 


is achievable using the codes L n *. both with ML decoding 
and naive lattice decoding. 

Proof: From equation (fl9b . we have that the hypotheses 
of Theorem 14.11 are satisfied (actually, only the weak law of 
large numbers was required). Consequently, any rate 


R < n log P + E h [log det H^H] — n log Cl + n log —- = 

/ p \ Th 

= log ( — j + E h [log det H^H] —n log Cl + n log = 


= E H 


p 

log det —H^H 
n 


- n log C L + nlog 


4 n 


is achievable. □ 

A similar corollary to Theorem 14.71 holds in the case n r < n: 


Corollary 5.6: Suppose that n r < n and that the fading 
process {Hi} is ergodic, stationary and isotropically invariant. 
Moreover, suppose that E [|logdetiTiTl|] < oo. Let L n ^ C 
AI n xnfc(C) be a family of 2n 2 fc-dimensional multiblock lattice 
codes such that det m m (L n> k) = 1 and Vol(L n ,fc)"^ ^ Cx. 
Then, any rate R lower than 


E 


H 


log det —7/7/1 
n 


— 2 n r —(n—n r ) log ■ 


n ire 

-hnlog-— 

n — n r nC l 


! We note that the capacity (per channel use) of the block fading MIMO 
channel of finite block length T with perfect channel state information at the 
receiver is independent of T (U, eq. (9)]. So the previous result still holds 
in the multiblock case. 
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is achievable using the codes L n .k both with ML decoding 
and naive lattice decoding. 

Remark 5.7: Using the same argument as in Remark 15721 
we can show that the achievable rate is within a constant 
gap from capacity, although this constant will depend on the 
fading model. 

Under the hypothesis Eh [Ilog det 72^721], we have 
|E h [log det 77^77] j < oo. The achievable rate is of 
the form R(P) = max (0, n log E + Eh [log det H^H] — c) 
for some constant c > 0. Let P m i n be the smallest value of 
P such that P(P) > 0 if P > P m i n . For P < P m i n we have 
that C(P) - R(P) = C(P) < C(P min ). 

Let Aj, * = 1,..., n be the (random) singular values of H. For 
P > Pmin, C(P) - R(P) = Er=l E H [log (l + ^t)] + C 
which is a strictly decreasing function of P and tends to c 
when P —>■ oo. This shows that the gap is uniformly bounded. 


VI. Achievable rates and error probability 

BOUNDS FOR I.I.D. RAYLEIGH FADING CHANNELS 


We now suppose that the coefficients of Hi are circular 
symmetric complex Gaussian with zero mean and unit variance 
per complex dimension, and that the fading blocks Hi are 
independent. In this case, the achievable rate can be computed 
explicitly, and we can prove that the error probability vanishes 
exponentially fast. 

Let ip(x) = ^lnr(x) denote the Digamma function. Then 
we have the following: 

Proposition 6.1: Let L rl k C M„ xn fc( C) be a fam¬ 
ily of 2n 2 fc-dimensional multiblock lattice codes such that 
det min(L n ,k) = 1 and Vol(L nifc )^ < C L . Then, over the 
( n,n r ,k ) multiblock channel, these codes achieve any rate 


R < Eh 


log det — 77^77 
n 


- n log C L + n log —, 
4 n 


where 


Eh 


log det —77^77 
n 


= log- 


P 


E V'W 


i=n r — 7J.+ 1 


( 20 ) 


Moreover, the error probability vanishes exponentially fast. 

Proof of Proposition 16.71 The first statement follows 
from Corollary 15.51 The next step is to prove equation (120b . It 
is well-known 14111421 1 that if H is an n r x n matrix with i.i.d. 
complex Gaussian entries having variance per real dimension 
1/2, the random variable det(//' //), corresponding to the 
determinant of the Wishart matrix II ' II , is distributed as the 
product 


F/,ro r — Z n „ — 


n+1 


Z n ~- 


n+2 * 


of n independent variables, such that Mj = n r — n+ 1,..., n r , 
2 Zj is a chi square random variable with 2 j degrees of 
freedom. The density of Zj is pzj{x) = X rj) -. We have 

i r°° 

EflnZj] = —— / x^~ 1 e~ x In a: dx = 

r(j) J o 

n r 

Mn t n r = E[ln Vn t n r ] = E f>{j) = E H [in det 77+] . 

j=n r —n+1 


Then if we consider the base 2 logarithm, we find 


Eh [log det 77+] = E[log V n , nr ] = 
_ Ei=„ r -„+i iPU) 


M, 


n,n r 


In 2 


= n log e n 


In 2 

A E ^n r -n +1 VU) 


which concludes the proof of equation ( l20t . 

In order to show that the error probability converges exponen¬ 
tially fast, by Remark 14.61 it is enough to show that we have 
exponential convergence in equation ( 1 1 4b . 

Consider a sequence of i.i.d. random variables In Vn ) lr , i = 
1,..., k, with the same distribution as lnV rainr . Using the 
Chernoff bound 114311 . given 6 > 0, Vi> > 0 we have 


P- 


M n n 1 X ' * TT ^ ^ I 

^logdet Hj Hi > — \ = 


In 2 


= P |M„ iTlr - - E In det H'\ H, > <5 j = 
= P j Af„,„ r - ^ E ln l+ r > <5 j < 

< e kHMn,n r S) (E[e- VlnVn ^]) k 


( 21 ) 


The tightest bound in (liHT l is obtained for vg such that 

E[— In V nt n r e~ vs ln W,n r ] = (S-M n , nr ) E[e“^ lny "+. 

Observe that 

Elz P = W)L 


r °° xi- x - v e~ x dx = r(j V) 


nZT-lnZj] = 


1 


^ ~ m 

r (j - v ) 


m ’ 

V j-l-v e -x l nxdx = 


( 22 ) 


m 

Thus we find 


-4>ti - v). 


(23) 


E [t 


—v In V n . 


n 


r}=E[V~l]= II E [V] = 

j=n r —n -\-1 

r (j - v) 


j=n r -n-\-l 

E [— \nV n , nr e~ vlnV ' 

n r 

= E E 

j=n r -n -\-1 


m 5 

] =E[-V~llnV n , nr ] - 

n r 

-In Zi [] Z T 


l=n r -n -\-1 


= E [l[nzr})n-z- v \nz J \ = 

j=n r —n +1 \l^j 


e n 

j=n r —n +1 \l^j 


r(i - v) | ry , 

T(irHvr*“" )= 


n r T'/7 A nr 

n 4 ^ e *u-v) 


T(l) 

l=n r — n+1 j=n r — n+1 
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Consequently, the tightest bound in (OH is achieved for vg 
such that 


n r 

5= Y, (m-W-vs)). (24) 

l=n r -n -\-1 


Note that as <5 —► 0, v's —t 0. The right-hand side in equation 
ED for v = Vs can be rewritten as 


e kv s {-5+Y^L nr _ n+1 4>(j)) 



r (i-vs)\ 

m ) 


k 


— e k(-vsS+J2”’Z nr -n+i( v s' , l } U)+ lnr U- v s)-inr(j))) _ 


— e k'Ej= nr .- n+ i(vs / <l>(j-v s )-lnr(j)+lnr(j-Vs)) 


using El. 

Recall that T(a;) is monotone decreasing for 0 < x < «o = 
1.461632... and monotone increasing for x > ao- Using the 
mean value theorem for the function In 1'() in the interval 
[i — vg, /] we get that for i = 1, vgtjt{l — t) 5 ) + lnr(l — vg) < 0, 
and for i > 2, vgipif — vs) < InT(i) — InT(z — vg). Thus, the 
exponent is negative both for n = n r and for n > n r . We can 
conclude that 


M, 


»=1 ) 


— kK n 


for some positive constant I\ n , nr: s. Since the bound holds 
V<5 > 0, using Remark l4~6l with p = , we find that the 

error probability tends to 0 exponentially fast for any rate 

R < n (log —e" 52i="r-n+i — log C L + log . □ 

V n An 


Corollary 6.2: Over the (n, n, k) multiblock channel, reli¬ 
able communication is guaranteed when k —> oo for rates 

R < n (^log — e" + log — log23 I0 ( 1_ "^G N ) 

\ n 2 n ) 

when using the multiblock code construction in Proposition 

El 


VII. Existence of asymptotically good lattices 


All of our capacity results depend on the existence of lattices 
with asymptotically good normalized minimum determinants, 
which was claimed in Section IIH-BI In this section we will 
prove this result. 

We will first recall the construction of single-block space- 
time codes from cyclic division algebras (see for example 
119]). Due to space constraints, we refer the reader to [^4] 
for algebraic definitions. 

Definition 7.1: Let K be an algebraic number field of 
degree to and assume that E/K is a cyclic Galois extension 
of degree n with Galois group Gal(E/K) = (o). We can 
define an associative /('-algebra 


A = {E/K, a, 7) = E ® uE ® u 2 E < 


®u n ~ 1 E, 


where u £ A is an auxiliary generating element subject to the 
relations xu = ua(x) for all x € E and «" = 7 e K*. We 
call the resulting algebra a cyclic algebra. 


Here K is the center of the algebra A. 

Definition 7.2: We call \J[A : K] the degree of the algebra 
A. It is easily verified that the degree of A is equal to n. 

We consider A as a right vector space over E and note that 
every element a = x 0 + ux\ + • • • + u n ~ l x n ^ 1 € A has the 
following representation as a matrix: 



( xo 

yo{x n -i) 

ya 2 (x n - 2 ) ■ 

• 7o-"“ 1 (:ri)\ 


Xi 

o(x 0 ) 

l<x 2 {x n -i) 

ja n ~ 1 {x 2 ) 

</>(a) = 

X2 

o(xi) 

(7 2 {xq) 

'ya n - 1 (x 3 ) 


V^n-l 

e{x n -2) 

a 2 (x n - 3 ) • 

■ (r n ~Hxo ) / 


The mapping 0 is called the left regular representation of 
A and allows us to embed any cyclic algebra into M n { C). 
Under such an embedding <j>{A) forms an mn 2 -dimensional 
Q-vector space. 

We are particularly interested in algebras A for which </>{a) 
is invertible for all non-zero a € A. 

Definition 7.3: A cyclic /('-algebra D is a division algebra 
if every non-zero element of D is invertible. 

If we assume that D is a division algebra, then cf> is an 
injective mapping to M r: (C) and every non-zero element in 
(j>{V) is invertible. However, 4>(D) is not a lattice. Therefore 
we will instead consider a suitable subset of V. 

Definition 7.4: A Z -order A in V is a subring of V having 
the same identity element as D, and such that A is a finitely 
generated module over Z which generates D as a linear space 
over Q. 

With the previous definition, the set <j>(A) is a matrix lattice 
that can be used for coding over a single space-time block. 

A generalization of the embedding <j> to the multiblock case 
was proposed in G2 fl3ll for division algebras whose center 
K contains an imaginary quadratic field. In this paper we 
consider a more general multiblock construction developed in 
m, which applies to any totally complex center K. 

We say that a degree 2k number field I\ is totally complex if 
for every Q-embedding fii :/(■—( C the image afiK) includes 
complex elements. The field K has 2k distinct Q-embeddings 
fit : K C. As we assumed that K is totally complex, each 
of these embeddings is part of a complex conjugate pair. We 
will denote by /3* the embedding given by x i-A fii{x). 

For each we can find an embedding ^ £ 4 C 

such that cii \k = fit- This choice can be made in such a 
way that ai\ k = A- We will suppose that the embeddings 
{ai,... ,ct 2 k} have been ordered in such a way that a; = 
a l+rn , for 0 < i < m. Let a be an element of V and A = 4>{a). 
Consider the mapping ip : A 1 —> M nxn k(C) given by 


a 1-4 {ai(A),... ,a k (A)), (25) 


where each ai is extended to an embedding ct; : M n (E ) °4 

M n (C). 

The following result was proven in [17, Proposition 5]: 

Proposition 7.5: Let A be a Z-order in D and ip the previ¬ 
ously defined embedding. Then 79 (A) is a 2fcn 2 -dimensional 
lattice in M nxnk { C) which satisfies 

det min ( 79 (A)) = 1, Vol(<?(A)) = 2 -*"V|d(A/Z)| 















11 


and 


%( A )) 


2 2 kn 2 \ 1/4 " 

RamJ 


Lemma 7.8: Every number field Kk in the Martinet family 
has ideals P\ and Pi such that 

N k/ q{Pi) < 23 fe/10 and N K/Q (P 2 ) < 23 fe / 10 . 


Here d( A/Z) is the Z-discriminant of the order A. It is a 
non-zero integer we can associate to any Z-order of V. We 
refer the reader to [44] for the relevant definitions. 

We can now see that in order to maximize the minimum 
determinant of a multi-block code, we have to minimize the 
Z-discriminant of the corresponding Z-order A. 

The first step to attack this question is to assume that A 
has some extra structure. Let Ok be the ring of algebraic 
integers of K. If we assume that A is also an Ok module, 
then the Ok -discriminant of A is well-defined 144], and will 
be denoted by d(A/0^). The following formula holds: 


d( A/Z) = N K/ q{d(A/0 K )){d K ) r 


(26) 


where N k /q is the algebraic norm in K and (Ik is the 
discriminant of the field K. 

In the case of fixed center K, d addressed the prob¬ 
lem of finding the division algebras with the smallest Ok~ 
discriminant, yielding the densest MIMO lattices. The main 
construction is based on the following result (Theorem 6.14 
in JH): 

Theorem 7.6: Let K be a number field of degree 2k and Pi 
and Pi be two prime ideals of K. Then there exists a degree 
n division algebra T> having an Ok -order A with discriminant 

d( A/Z) = (N K/Q (P 1 )N K/Q (I* l )) n t n - 1 '>(d K ) n \ (27) 

Theorem 17.61 suggests that in order to build families of 
(n, n, k) multiblock codes with the largest normalized min¬ 
imum determinant, we should proceed in two steps: 

a) choose a sequence of center fields K of degree 2k such 
that their discriminants dK grow as slowly as possible; 

b) given the center K, choose an algebra T> satisfying (l27l >. 
where Pi and P 2 are the prime ideals in K with the 
smallest normal 

We now discuss the choice of a suitable sequence of center 
fields. The following theorem by Martinet [45] proves the 
existence of infinite sequences of totally complex number 
fields K with small discriminants d k . As we will see in the 
following, choosing such a field as the center of the algebra 
D is a key element to obtain a good normalized minimum 
determinant. 

Theorem 7.7 (Martinet): There exists an infinite tower of 
totally complex number fields {Kk} of degree 2k, where 2k = 
5 • 2 2+t , such that 

I d Kk \^=G, (28) 


for G w 92.368. 

The following Lemma shows that the number fields in the 
Martinet family have suitable primes of small norm yielding 
a good bound in Theorem 17.61 

2 However, we note fa that a priori there may be a tr ade-off between 
these two choices, so that minimizing the two terms in CD separately may 
be suboptimal. 


Proof: Every field Kk has a subfield F = 
0(cos(27r/ll), s/2, \/ —23), where [P : Q] = 20 (see 
for example 0461 p. 395]). The field F has prime ideals B\ 
and Bi such that AA/q(P,:) = 23. Let us now suppose that 
Pi and P 2 are such prime ideals of Kk that P, n O p = II,. 
Transitivity of the norm then gives us that 

N Kk/Q (Pi) = Afp/q(N Kk /F{Pi)) < Afp/q (P.;) 2fc / 20 . □ 


Armed with this observation, we can finally prove Proposition 

[L8l 

Proof of Proposition I3.<S1 

Suppose that we have a degree 2k field extension K in 
the Martinet family of totally complex fields such that (l28l) 
holds. We know that this field K has some primes P\ and 
Pi such that Nk/qSPi) < 23 fc / 10 and N K/Q {P 2 ) < 23 fc / 10 . 
Then, there exists a central division algebra T> of degree n 
over K, and a maximal order A of T>, such that 

d( A/Z) = (N k/Q (P 1 )N K/Q (p i ))« n - 1 '>(d K ) nl ‘ < 

< (23 fc / 5 )( ri ( ri — x ))(C7 2fc ) ri2 □ 


Remark 7.9: We note that the number field towers in The¬ 
orem 17.71 are not the best known possible. It was shown in 
0 that one can construct a family of totally complex fields 
such that G < 82.2, but this choice would add some notational 
complications. 

VIII. Corollaries to single antenna fading 

CHANNEL 

The single antenna fast fading channel is one of the special 
cases of the general channel model ©. It is particularly 
illuminating as the connection to the classical AWGN lattice 
coding is most striking. In this case the abstract matrix lattices 
of Section IVIII correspond to simple number field codes that 
have been studied for twenty years 0. Due to the familiarity 
and simplicity of this model we can most easily compare our 
work to previous research on the topic. 

In the single antenna case the channel model © gets 
simplified to 

Di = hi- Xi+Wi , (29) 

where Xi are the transmitted symbols, and V* = 1 ,... ,k, Wi 
are i.i.d. complex Gaussian random variables with variance 
a 2 = a 2 = \ per real dimension and { h ,} is some complex 
fading process such that JA =1 ^ log \hi\ 2 converges in prob¬ 
ability to some constant when the number of blocks k tends 
to infinity. 

This scenario has received considerable interest in the case 
of an i.i.d. complex Gaussian fading process {hi }, and several 
works have focused on the design of lattice codes for this 
model 0 0. The analysis of the union bound for the 
pairwise error probability for a lattice code L C C k leads to a 
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design criterion based on the maximization of the normalized 
product distance 


Nd B 


i{L) = inf 


nit \xi 


xeL\{o} Vol(L )2 


Note that the normalized product distance is a special case (for 
n = 1) of the normalized minimum determinant in <[6]». Most of 
the works in the literature have focused on the optimization 
of the product distance for lattice signal constellations with 
a fixed number of blocks k\ few authors M l49h have also 
studied the upper and lower bounds for Nd ?vmm over all 
lattices when k grows to infinity. 

However, there has been no general consensus on whether 
significant gain could be achieved from coding over an exten¬ 
sive number of fading realizations. For example the authors in 
[50] state that: “increasing the diversity does not necessarily 
increase to the same extent the performance: in fact, the 
minimum product distance decreases and the product kissing 
number increases. Simulations show that most of the gain is 
obtained for diversity orders up to 16”. In fact, the analysis 
of the distribution of pairwise errors in the union bound as in 
El shows that the product kissing number lEsIl . or number 
of worst case occurrences, will grow fast and a priori might 
eat away the product distance gain. However, this issue seems 
to be due to the suboptimality of the union bound rather than 
to the codes themselves. 

In fact, let us consider an infinite family of 2/c-dimensional 
lattices L k C C fc with normalized product distance satisfying 
(Nd p m i n (Lfc)) 2 / fe > c, for some positive constant c. 

According to Theorem 14.11 and Remark 14.31 we then have 
the following. 

Corollary 8.1: Any rate It 


R < E/i 


\og 2 P\h\ 2 +log : 


7re 


log 2 


is achievable with the family L k of lattices over the fading 
channel ( [29b . 

This result proves that indeed we gain by coding over an 
increasing number of blocks, assuming that we have a family 
of lattices Lk with the described product distances. According 
to Proposition 13.61 the condition (Nd Pim i n (£fe)) 2 ^ fe > c 
implies that l'ho(Tfc) > fee. It reveals that families of lattice 
codes with large product distance do not only have large 
Hermite invariants, but also that the Hermite invariants of the 
faded lattices are as large as well. Thus, the product distance 
is not only relevant in capacity considerations or in the high 
SNR scenario, but also plays a role when coding over a finite 
number of fading realizations for low SNR. 


A. Approaching capacity with number field codes 

Using the normalized product distance as a code design 
criterion led to lattice constructions based on number fields 
in [52, [221 48, 18]. However, none of these works considered 
capacity questions. 

Let us now show how the construction in Proposition 17.51 
when specialized to the single antenna case, is just the standard 
method used to build lattice codes from number fields 


and how this method can be used to approach capacity in fast 
fading channels. 

Let K/Q be a totally complex extension of degree 2k and 
{or,..., erfc} be a set of Q-embeddings, such that we have 
chosen one from each complex conjugate pair. Then we can 
define a relative canonical embedding of K into C n by 

ip(x) = (<7i (x),.. ,,a k (x)). 

The ring of algebraic integers Ok has a Z-basis W = 
{wi,...,W 2 k} and p(W) is a Z-basis for the full lattice 
<p(0 K ) in <C k . 

Proposition 17.51 now simplifies to the following. 

Corollary 8.2: Let ip be the previously defined embedding 
and K a degree 2k totally complex number field. Then 1 p(Ok) 
is a 2fc-dimensional lattice in C k which satisfies 

det m i n (^(Ok)) = 1, Vol(<y5(C7ft-)) = 2~ k y/\dx\ 


and 


S(v(Ok)) 


-L 


2k \ V 4 


\l dic\ 


Using Martinet’s family of fields K k from Theorem [777] and 
setting Li^k = p{OK k ) we have 

Vol(Li |fc ) < and det min (Li, fc ) = 1, 

where G ~ 92.368. Specializing to the case where the fading 
process is i.i.d complex Gaussian we have that any rate 

R < log 2 (Pe -7 ) - log 2 ■ (30) 

where e -7 = Eh[ log |/t,;| 2 ], is achievable. 


B. Known bounds on discriminants and Hermite invariants 
Equation ( l30l > reveals that the codes based on the Martinet 
family have a rather large gap to capacity. However, the right- 
hand side of (f30b is just a lower bound on the maximum 
achievable rate with lattice codes, and might be improved with 
a better error probability estimate and/or a better choice of the 
lattice sequence. 

An upper bound for the maximum achievable rate using this 
approach can be derived from a lower bound for the discrim¬ 
inant. The Odlyzko bound [53] states that when k —> 00 we 
have that 1 1 / 2fe > 22.3. If it were possible to reach this 
lower bound with an ensemble of lattice codes, then any rate 
R satisfying 

/ 44.6\ 

R < log 2 (Pe 7 ) - log 2 ( — ) , (31) 


would be achievable. For small values of k, there exist number 
fields having considerably smaller root discriminants. Table |T| 
El lists the best known root discriminants for totally complex 
number fields of degree 2k. The first four values are known 
to be optimal. 

We note that even equation m does not represent an 
absolute limit for the rates that are achievable with lattice 
codes and not even with algebraic lattices arising from number 
fields, and does not mean that the performance of algebraic 
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TABLE I 

Best known root discriminants for totally complex number 

FIELDS K OF SMALL DEGREE 2 k . 


k 

\d K \^ 2k 

1 

1.732.. 

2 

3.289.. 

3 

4.622.. 

4 

5.787.. 

5 

6.793.. 


codes will always be bounded away from capacity. In fact, 
as seen in Corollary 18.11 we are only interested in the nor¬ 
malized product distance of the lattices under consideration. 
For example, instead of considering the image of the ring of 
integers Ok under the embedding <p, one can use an ideal of 


this ring of integers (see 154 1491 Li!]) or more generally any 


lattices with good normalized product distance. 

The Minkowski-Hlawka theorem provides a non¬ 
constructive proof of the existence of 2fc-dimensional 
lattices Lk C C k having Hermite invariants h (Lk) ~ — [6]. 
If it were possible to obtain also hc(Lk) ~ or equivalently 
(Ndp, min (Lfc)) - (^) fc/2 , then all rates satisfying 


R<Eh 


\og 2 p\hy 


. 4 , 1 

- log 2 -b log2 — 

7re 7re 


= E h \og 2 P\h\ 2 — 2 


(32) 


would be achievable with this family of lattices. However, 
we do not know if this is possible. The two bit gap to 
Eh [ log 2 P \h\ 2 ] in (l32l > would be exactly the same that is 
obtained in the AWGN case when using the hard sphere 
packing approach Chapter 3]. Still this two bit gap is not a 
fundamental limit of the performance of lattice codes but likely 
an artifact of the suboptimal method to analyze the error. 


IX. Geometry of numbers for fading channels 

In the previous sections, we have shown that the normalized 
minimum determinant provides a design criterion to build 
capacity-approaching lattice codes for block fading multiple 
antenna channels. Let us now see how this approach fits into 
a more general context and can be regarded as a natural 
generalization of the classical theory of lattices for Gaussian 
channels. Finally we show how the code design problems, both 
in Gaussian and fading channels, can be seen as instances of 
the same problem in the mathematical theory of geometry of 
numbers. 

Consider a lattice L C C k having fundamental parallelotope 
of volume one and define a function /i : C k —> R by 

fi(xi, ...,Xk) = \xi\' 2 + \x 2 \ 2 H-F \x k \ 2 - (33) 

The real number h(L) = inf^gj^ x ^o f\ (x) is then the Hermite 
invariant of the lattice L. Let us now denote with Ck the set 
of all 2fc-dimensional lattices in C k with volume one. 

Suppose that we have an infinite family of lattices Lk G Ck 
with Hermite invariants satisfying > c, for some positive 


constant c. As stated in the beginning of Section [TV] then all 
rates satisfying 


R < log 2 (P) - log 2 



+ log 2 c, 


are achievable in the complex Gaussian channel, with this 
family of lattices. The Hermite invariant h {Lk) now roughly 
describes the performance of the lattice Lk and can be used 
to estimate how close to the capacity a family of lattices can 
get. This relation is one of the key connections between the 
theory of lattices and information theory and has sparked 
a remarkable amount of research. 


Let us now see how our results can be seen as natural 
generalizations of the relation between Hermite invariant and 
capacity. 

Let us consider 2/.:-dimensional lattices L C C k in Ck and 
the form 


f 2 (x 1 ,X2,..-,Xk) = |anx2 ■ ■ - Xk\- (34) 


Then Nd p ni i n (L) = f 2 (x), is the normalized 

product distance of the lattice L. 

Assume that we have an infinite family of lattices 
Lk G Ck with normalized product distance satisfying 
(Nd Pim in(L ra )) 2 / fc > c, for some positive constant c. As seen 
before we have that all rates satisfying 


R<E h 


\o g2 p\hy 


- log 2 -b log 2 c 

7re 


are accessible with this family of lattices with zero error 
probability over the Rayleigh fast fading channel. 

We denote with C^„ k) the set of all 2n 2 fc-dimensional 
lattices in the space M nxn k(C) with volume one. Given 
(Xi,X 2 , ..., Xk) G M nX kn{ C), we consider the function 


k 

f 3 (X 1 ,X 2 ,...,X k )=l[ |detpQ)|. 

i=1 


For L G £( n ,fc)> we have 


S(L) = 


inf 

X£L k ,Xjt 0 


h(X). 


If Lk C Af nX fc ra (C) is a family of lattices with the property 
that 5(Lk) 2 / kn > c then according to Remark [4.31 any rate 
satisfying 


R < Eh 


p 

log 2 det —H^H 
n 


+ n log 2 — + n log 2 c 


is achievable with the lattices Lk- 

We can now see that the normalized minimum determinant 
and product distance can be regarded as generalizations of 
the Hermite invariant which characterize the gap to capacity 
achievable with a certain family of lattice codes. 

A natural question is how close to capacity we can get with 
these methods by taking the best possible lattice sequences. 

The Hermite constant H(k) can now be defined as 


fT(2fc) = sup{h(L) | Lg £ (1 , fc) }. (35) 

In the same manner we can define 


Ndp, min (fc) = sup{Nd p , min (L) | L G £(i, fc )}. (36) 
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and 

6{k,n) = sup{S(L) | L e £(„,&)}■ (37) 


Each of these constants now represents how close to capacity 
our methods can take us. Any asymptotic lower bound with 
respect to k will immediately provide a lower bound for the 
achievable rate. Just as well upper bounds will give upper 
bounds for the rates that are approachable with this method. 


The characterizations of achievable rates using lattice codes 
have now been transformed into purely geometrical questions 
about the existence of lattices with certain properties. The 
value of the Hermite constant H(k), for different values of 
k, has been studied in mathematics for hundreds of years and 
there exists an extensive literature on the topic. In particular 
good upper and lower bounds are available and it has been 
proven that we can get quite close to Gaussian capacity with 
this approach [6, Chapter 3]. 

In the case of the product distance, this problem has been 
considered in the context of algebraic number fields and some 
upper bounds have been provided. As far as we know the best 
lower bounds come from the existence results provided by 
number field constructions PiJ and fill - 

The properties of S(k, n) have been far less researched in the 
literature. Simple upper bounds can be derived from bounds 
for Hermite constants as pointed out in [33] and lower bounds 
are obtained from division algebra constructions as described 
in this paper, but the mathematical literature doesn’t seem to 
offer any ready-made results for this problem. 

However, all three questions can be seen as special cases of the 
problem of finding the minimum of a homogeneous form over 
a lattice in the mathematical theory of geometry of numbers 
[35]. Let us now elaborate on the topic. 

Definition 9.1: A continuous function F: M nX kn(C) —> R 
is called a homogeneous form of degree a > 0 if it satisfies 
the relation 


\F(aX)\ = \a\°\F(X)\ (Va G R,VX G M„ xfen (C)). 

Let us consider the body S(F) = {X | X e 
M nxkn {C),\F{X)\ < 1}, and a 2 kn 2 dimensional lattice L 
with a fundamental parallelotope of volume one. 

We then define the homogeneous minima \(F. L) of F with 
respect to the lattice L by 

A (F, L) = (inf{A| A > 0, dim(R(AS'(F) n L)) > 1}) CT , 

where M.(XS(F) IT L) is the R-linear space generated by the 
elements in A S(F) IT L. This allows us to define the absolute 
homogeneous minimum 

X(F) = sup A (F,L). 

Vol(L) —1 

We can now see that all of our forms /i, fi and are 
homogeneous forms. For the Hermite invariant we have a = 
2 , for the product distance a = n, and for the normalized 
minimum determinant a = n 2 k. We can also easily see that 
the constants (l35l l. d36l > and d37l > are absolute homogeneous 
minima of the corresponding forms. 

These results suggest that there is a very general connection 
between information theory and geometry of numbers for 


different channel models. It seems that given a fading channel 
model, there exists a form whose absolute homogeneous 
minima provide a lower bound for the achievable rate using 
lattice codes. 

Remark 9.2: The definitions for the geometry of numbers 
given in this section were stated for lattices in the space 
ALraxnfc(C), while usually the definitions are given in the space 
R m . This is however, just to keep our notation simple. The 
space M nxn k(C) can be identified with the space R 2 ™ k and 
we could have given the definitions also in the traditional form 
using this identification. 

X. Discussion and questions for further research 

In this work we proved the existence of lattice codes achiev¬ 
ing constant gap to capacity in ergodic fading channels. Unlike 
the case of existence results based on random coding, our 
finite codes are always built from the same family of lattices, 
irrespective of the SNR and even of the fading statistics. 
Hence, using the minimum determinant as a design principle 
leads to extremely robust codes. In particular division algebra 
and number field codes have this robustness property. 

However, our codes still have a considerable gap to capacity 
and further research is needed. Let us now point out a few 
directions this research can take next. 

In the case of single user channels the clearest goal is 
to improve our methods and close the gap to capacity. We 
note that this gap depends on several factors. First of all, the 
normalized minimum determinant affects the value of the gap. 
Second, our bound for the error probability is based on sphere 
packing and thus is suboptimal. 

Thus, the possible improvements to our construction are 
two-fold. In the first place, one could try to find families of 
lattices L n C M n ra fc(C) with larger normalized minimum 
determinant, for instance by replacing the centers in our 
constructions with families of number fields having smaller 
discriminants. One can also consider more general examples of 
lattices than those arising from orders in division algebras: for 
example ideals of orders, or in the case of number field codes, 
ideals of the ring of algebraic integers. In the second place, 
in this paper we have not considered the issue of shaping. 
Improving the shaping properties of our lattices might lead to 
a better error probability bound. 

Another approach is to relax our minimum determinant code 
design criterion. Our codes are extremely robust and quite 
universal in the sense that they respond very well to any non 
pathological fading realization. This universality is of course a 
strength, but it could also lead to a situation where the codes 
are rather good for every channel, but not optimal for any. 
If we fix a channel model, it may be possible to weaken 
the design principle. This might allow us to consider larger 
ensembles of lattices and possibly to close the gap to capacity 
in this fixed channel model. 

Our work was about explicit code constructions in the spirit 
of classical sphere packings [@]. However, it seems that even 
the existence of capacity achieving lattice codes in fading 
channels is an open question (see Section Il-Ab . In the case of 
AWGN channel this question was solved only quite recently in 
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[| 25ll by assuming that the receiver and transmitter have access 
to a common source of randomness. One of the key elements 
in the achievability result of i25ll is the Minkowski-Hlawka 
theorem, that can be used to prove the existence of lattices with 
certain properties. It is then a natural idea to generalize this 
approach by proving an analogue of the Minkowski-Hlawka 
theorem for fading channels. It seems to us that this problem 
is non-trivial. 

In this paper we have considered block fading MIMO 
channels, but we hope that the methods developed here can 
be applied also in a more general setting. Let us now sketch 
an outline for possible generalizations. 

The reduced Hermite invariant is a natural analogue of the 
classical Hermite invariant for fading channels. This concept 
can likely be generalized to other fading channel models, such 
as for example intersymbol interference channels. Given a 
fading channel we can ask what would be the group (or set) 
G that would represent the action of the channel, and define 
the corresponding reduced Hermite invariant Iiq. The next 
question is then to find lattices that would maximize this value. 
In the case of the block fading channel, the problem was made 
more accessible by Proposition [376] where we proved that ho 
can be seen as the minimum of a certain homogeneous form. 
This line of thought suggests a general approach to turn the 
chase for capacity into a problem in geometry of numbers for 
different channel models. It also raises several questions. For 
example we can ask which are the channel models where this 
approach can be applied and for which groups G the reduced 
Hermite invariant corresponds to some homogeneous form. 

Finally, the lattice codes proposed in this paper could have 
applications to other problems in information theory, such as 
coding for multiple access fading channels and for information 
theoretic security. 


Appendix 


A. Proof of Theorem \4.7\ 

With a similar approach as in the proof of Theorem 14. II we 
consider the following upper bound: 


p e <p<! \\w\\ 2 >[^f 


< 


< 


\\wtf 

knn r 


> 1 + e > +1 


,r di 

\Aknn, 


< 1 + e 


(38) 


The first term in equation (f38l> tends to zero exponentially fast 
when k — > 00 since 2||kF|| 2 ~ % 2 (2 knn r ). We now focus 
on the second term in equation (l38l >. and begin by finding 
a lower bound on the minimum distance du in the received 
constellation. 

For all i G {1,..., k}, let A ij, j = 1,..., n be the singular 
values of // ■' II, with 


0 — Aj ? i — * — A i : n—n r ^ f f A i : m 

and lij the singular values of XiXJ with 


h, 1 > h ,2 > • • • > h,-, 


Using the mismatched eigenvalue bound [55j,[20[], we have 

n n 

IIMill — 'y = 

j —1 j—n-n r +1 

Consequently, we find that 


d 2 H > a 2 _ min _ Ajkj > 


xeL„, fc \{ 0} 


2=1 j=n— n r +l 


^ OL n r k | j | (A i^j lij ) nrk 

2=1 j=n—n r -\- 1 

Using the NVD property of the code, we get 

k n k 

nn^niMi 2 ^ 

2=1 j = l 2=1 

Therefore, we have the lower bound 

k n / k n—n r \ — * 

n n a > (n n a j 

\i=i 1 / 

k n-nr \ -(n-n r )k 


(39) 


> 


2=1 j=n— n r + l 


> 


> 


1 


(n - n r )k A A 

^ 2=1 7=1 


li d 


> 


Pn 


2 (n — n r ) 


-k(n—n r ) 


where we have used the arithmetic-geometric mean inequality 
and the power constraint a 2 ||X|| 2 = a 2 l-‘,j ^ Pkn. 

Replacing the previous expression in ( l39l ). we obtain 


d 2 H > 


a 2 n r k Till n"=„-„ r +i K, : 


( Pn V 

l a 2 (n—n r ) J 


= (a 2 f 


n — n r 
Pn 


,/e n det ( Hj H, 




i=i 


The second term in (1381) can thus be upper bounded by 


i=l 


( n ) 

I nr 1 

( p 

\a 2 ) 


\n — n r J J 


■ ^ log det HiHj < log 


t ^^( 4(1 + e))"-n"P" 


2 = 1 


a 2n (n — n r ) r 


By hypothesis the weak law of large numbers (IT5l i holds, 
i.e. i Y^i=i log det HiHj —> /i as k —> 00 . Thus, the error 
probability will vanish provided that for sufficiently large k, 

, (4(1 + e ))nr n n pn-n r 

lQ g -2 ~l -- < d 

a zn (n~n r ) n Ur 


Q n 2 k p 

Recalling that a 2 > -ff —, the condition can be rewritten as 

2 ~C L 

R < [i + n r logP — n r log 4(1 + e) — n log uCl+ 

, logCn.fe 


nk 


+ (n — n r ) log(n — n r ) 
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Using Stirling’s approximation (fT3] >. for large k we have 

log C’n k , , 1 , r, 2i 

- - - ft! n log 7 re — n log n — —— log 27 m k 

nk 2 nk 

Asymptotically, we find that any rate 
P . n 2 C L 


R < ii + n r log 


4(1 + 6) 


— n log ■ 


7 re 


+ (n — n r ) log (n — n r ) 


is achievable. Since this is true for all e > 0, this concludes 
the proof. □ 
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